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1.  Introduction 

In  the  Al  literature,  the  topics  of  qualitative  reasoning  and  causality  have  been  closely 
interwoven,  and  have  sometimes  been  treated  as  though  there  was  an  essential  connection 
between  them.  In  fact  however,  the  two  topics  are  quite  separable.  Qualitative  reasoning  is 
mainly  used  to  solve  problems  in  comparative  statics:  if  there  is  a  change  in  the  value  of  an 
exogenous  (independent)  variable  of  a  system  of  variables,  what  changes  will  take  place  in  the 
endogenous  (dependent)  variables  when  the  system  reaches  its  new  equilibrium?  Causal 
reasoning,  on  the  other  hand,  is  mainly  used  to  explain  how  a  system  works:  the  mechanisms  of 
causal  connection  among  its  parts  and  the  way  in  which  these  causal  connections  propagate 
effects  through  the  system. 

Of  course  the  two  methods  of  analysis  can  be  combined.  The  causal  ordering  of  the 
variables  of  a  system,  together  with  the  algebraic  signs  associated  with  the  causal  links  provide 
sufficient  information  to  carry  out  some  (but  not  ail)  qualitative  analyses.  However,  given  any  set 
of  equations,  causal  or  not,  that  describe  a  system,  and  given  the  signs  of  the  coefficients  in 
these  equations,  qualitative  analysis  can  be  carried  out  without  concern  for  the  causal  ordering  of 
the  variables. 

In  previous  papers,  two  of  us  have  explicated  the  formal  foundations  of  causal  ordering  for 
dynamic  systems,  static  systems,  and  systems  containing  a  mixture  of  dynamic  and  static 
mechanisms.  We  have  also  treated  the  problem  of  aggregation  and  disaggregation  of  causal 
systems,  and  have  shown  both  how  to  carry  out  these  processes  and  what  their  consequences 
are  for  the  descriptions  of  the  systems  [4, 3, 5]. 

In  this  paper,  we  propose  to  provide  comparable  foundations  for  qualitative  analysis,  and  in 
particular,  to  show  how  qualitative  reasoning  arises  naturally  and  simply  out  of  the  structure  of 
systems  of  algebraic  equations  and  ordinary  differential  equations.  Our  goal  is  to  sum  up  and 
explicate  the  practices  of  researchers  in  many  fields  who  use  qualitative  reasoning,  and  thereby 
to  gain  an  understanding  of  the  formal  assumptions  and  mechanisms  that  underlie  this  kind  of 
analysis.  There  is  no  new  mathematics  in  what  we  have  to  say.  Our  reason  for  saying  it  is  that 
there  does  not  appear  to  exist  any  readily  available  general  and  systematic  account  in  the  Al 
literature  of  the  simple  mathematics  that  is  implicit  in  what  people  are  doing  when  they  carry  out 
qualitative  analyses  of  systems  of  interconnected  variables.1  The  extensive  discussions  of 

'Perhaps  the  bast  treatments  of  the  subject  -  al  now  very  old  --  are  to  be  found  in  Lotka  [6],  Evans  [2],  and  Samuel  son 
(8).  the  first  of  which  applies  the  analysis  mainly  to  biology,  the  last  two  to  economics. 
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qualitative  reasoning  to  be  found  in  the  Ai  literature  are  not  rooted  in  the  familiar  formalisms  of  the 
algebra  and  calculus  of  real  numbers  but  postulate  new,  independent  formalisms. 

Since  it  has  not  proved  easy  to  make  these  new  formalisms  consistent  and  rigorous,  it 
would  be  very  convenient  if  we  could  simply  take  advantage,  for  this  purpose,  of  familiar 
mathematical  systems  that  are  already  in  place.  We  will  show  that  this  can  indeed  be  done,  and 
that  standard  formalisms  can  be  used  to  explicate  the  practices  of  qualitative  reasoning. 

2.  Qualitative  Reasoning  Based  on  Ordinal  Relations 

Qualitative  reasoning  treats  quantities  that  are  scaled  only  ordinally,  hence  are  defined  only 
up  to  arbitrary  strictly  monotonic  transformations.  Greater  and  less,  warmer  and  cooler,  taller  and 
shorter  are  ordinal  terms  that  can  be  used  to  describe  the  qualitative  relations  between  pairs  of 
values  of  the  corresponding  ordinal  variables  (amount  temperature,  height).  On  the  other  hand, 
terms  like  "twice,"  "half,"  "four  more  than"  imply  cardinality  and  cannot  be  applied  to  ordinal 
variables.  If  ardor  is  an  ordinal  variable,  then  one  lover  cannot  be  twice  as  ardent  or  half  as 
ardent  as  another  -  only  more  or  less  ardent. 

Of  course,  the  fact  that  cardinal  properties  cannot  be  attributed  to  ordinal  variables  does 
not  prevent  us  from  attributing  ordinal  properties  to  cardinal  variables.  Even  in  the  presence  of  a 
Celsius  thermometer,  one  day  can  simply  be  warmer  than  another,  or  cooler.  In  particular,  it  may 
often  be  convenient  to  adopt  a  particular  cardinal  scale  to  represent  an  ordinal  variable,  restricting 
ourselves  to  just  the  ordinal  properties  of  the  scale  and  renouncing  any  interpretation  of  its 
cardinal  properties.  Thus,  in  the  (ordinal)  Mohs  scale  of  hardness,  cardinal  numbers  are 
assigned:  for  example,  2  for  gypsum,  4  for  fluorite.  One  can  say  that  fluorite  is  harder  than 
gypsum,  but  not  that  it  is  twice  as  hard.  Still,  it  may  be  convenient  to  employ  the  numerical  scale, 
provided  that  we  remember  to  attach  meaning  only  to  the  ordinal,  and  not  to  the  cardinal 
properties  of  the  scale. 

Ordinal  relations  are  invariant  under  arbitrary  strictly  positive  monotonic  transformations  of 
the  scale  of  measurement,  while  cardinal  relations  are  invariant  only  under  affine  transformations. 
Although  ordinal  relations  are  weaker  than  cardinal  relations  (the  former  are  implied  by,  but  do 
not  imply,  the  latter),  they  have  the  advantage  that  they  can  often  be  established  empirically  when 
we  do  not  have  enough  information  to  establish  values  on  a  cardinal  scale.  Thus,  when  I  step  out 
of  doors,  I  may  know  that  it  is  warmer  than  it  was  yesterday  without  knowing  the  exact  F <u  irenheit 
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or  Celsius  temperature  on  either  day. 

If  we  ari  considering  an  equation  connecting  two  variables,  y~ff,x),  we  may  know  that  as  x 
increases  y  win  increase  (dy/dx> 0)  without  knowing  the  exact  functional  form  of  the  relation,  much 
less  the  values  of  numerical  parameters.  Or,  even  if  we  can  know  the  form  of  the  function  and 
the  parameter  values  for  one  particular  situation,  we  may  be  interested  in  reasoning  about  the 
equation  generally,  over  a  whole  range  of  situations  of  which  this  one  is  only  a  special  case,  but 
in  all  of  which  y  remains  some  strictly  monotonically  increasing  function  of  x.  In  both  of  these 
cases,  we  will  want  to  use  the  methods  of  qualitative  (ordinal)  reasoning.  Let  us  now  state  these 
matters  a  little  more  formally. 

First,  a  few  words  about  notations:  we  will  sometimes  use/,  to  denote  the  derivative  of/ 
with  respect  to  x.  If/is  a  function  of  x  alone,/,  is  the  total  derivative  of/with  respect  to/.  When / 
is  a  function  of  several  variables  including  x,  /,  is  a  partial  derivative.  Whenever  this  notation  is 
used,  the  context  should  make  it  dear  which  is  the  case. 

2.1.  Invariance  under  Monotonlc  Transformation 

If  a  variable  is  meant  only  to  define  the  ordering  of  some  items  along  a  scale,  we  style  the 
variable  ordinal.  This  variable  can  be  replaced  by  any  other  that  does  not  disturb  the  ordering  by 
transposing  items.  Thus,  if  x(x>0)  is  a  real  variable,  but  we  are  interested  only  in  its  ordinal 
properties,  we  may  replace  it  by  x2  or  e*  or  logx,  without  changing  these  properties,  for  if  x,  >x; 
then  x2,.  >x2  ,  >  e*j  and  logxi  >  logxj. 

More  precisely,  the  ordering  of  a  set  of  elements  by  a  variable  associated  with  that  ordering 
is  invariant  under  any  positive  strictly  monotonic  transformation  of  the  variable.  Consider  xzX 
and  ye  Y.  Then  ] \x)-X-*Y,  is  a  positive  strict  monotonic  transformation  if  y=J{x)  >y;  */(xp 
whenever  x(  >xj  and  vice  versa.  If  a  set  of  objects  has  been  ordered  by  the  assignment  of  values 
of  x  then  replacement  of  all  these  values  by  the  corresponding  values  of  y  will  not  change  the 
ordering.2 

For  simplidty,  we  will  confine  ourselves  to  variables  that  are  defined  over  the  reals.  Let  x 
and  y  be  a  pair  of  such  variables,  and  let  us  suppose  that  they  are  positively  monotonically 
related.  That  is,  there  is  a  function /»,/*-►/?,  such  that  x,  >  xy  iff  y- */(x1)  >  y;*/(x;).  Now  we 

*ln  what  follow*,  wa  will  omit  the  adjective  ‘strict*  from  our  charactarization  of  monotonic  transformations,  but  for 
simplicity,  wa  will  limit  oursafvas  to  tha  rotation,  >,  of  strict  inequality  instead  of  using  the  weaker  relation,  >,  ‘greater  than 
or  equal  to.' 
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subject  ooth  x  and  y  to  positive  monotonic  transformations,  g(»)  and  h(»)  respectively,  so  that 
*•$(*)  and  wmh(y).  Since  g  and  h  are  monotonic  we  have  2t  >  z;  iff  xi  >  xJt  and  w(  >  w  iff  y{  >  yy 
But  since  x  and  y  are  positively  monotonically  related  xi  >  x;  iff  w(  >  wy  This  is  readily 
generalized  to: 

Theorem  1:  If  two  variables  x  and  y,  are  positively  (negatively)  monotonically 
related,  and  if  x  is  transformed  to  2  and  y  to  w  as  2 «  g(x)  and  w  «  h(y),  where  g  and  h 
are  arbitrary  positive  monotonic  transformations,  then  w  and  2  are  positively  (negatively) 
monotonically  related. 

Proof:  Since  g  is  a  monotonic  transformation 
X,  >  *;<-=>*,  >  *>• 

Similarly, 

y,  >  y/e— >w.  >  wj, 

since  h  is  monotonic  transformation.  Given, 
x,.  >  x<— >y,  >  yjt 

it  follows  that 

2-  >  Zj9  <**>x,  >  Xj,  -+y}  >  y^,  -*wi  >  Wj.  •  ■  ■  • 

For  example,  suppose  we  have  observed  empirically,  using  ordinal  measures,  that 
pressure  grows  with  temperature.  Then  if  we  replace  our  pressure  scale  with  a  new  ordinally 
equivalent  scale  (e.g.,  /**.  where  /*•  is  the  square  of  pressure),  and  our  temperature  scale  with  an 
ordinally  equivalent  scale  (e.g.,  7>,  where  T*  is  the  logarithm  of  temperature),  it  follows  that  /*• 
grows  with  T*. 

Moreover,  monotonic  relations  between  ordinal  variables  are  transitive  in  the  following 
sense: 

Theorem  2:  If  x  varies  positively  with  y  and  y  varies  positively  (negatively)  with  z, 
then  x  will  vary  positively  (negatively)  with  z.  The  latter  relation  will  be  invariant  under 
arbitrary  positive  (negative)  monotonic  transformations  of  x  and  2. 

Proof:  Given  y  >  0,  and  2  -  g(y),  gy  >  0, 

dx  dydx 

Therefore  monotonicity  is  transitive.  Note  this  holds  for  £ ,  $  <  0. 

dx  ay 

Let  v*A(x),  h *x(z)  such  that  vx,  wt  >  0. 

dw  _dw  dz  _dw  dz  dy  dx  ^  ^ 
dv  dz  dv  dz  dydx dv 

Note  that  £>0  since  —  >  0.  This  is  also  true  for  v„  w,  <0.  •  •  • 

d*  dx  x  1 
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3.  Continuous,  DlffSrsntiabls  Functions 

If  functions  of  ordinal  variables  can  be  represented  by  functions  of  cardinals  that  are 
continuous  and  possess  appropriate  derivatives,  then  for  monotonic  transformations  that  are  also 
continuous  and  possess  appropriate  derivatives,  we  can  apply  the  usual  operations  of  the 
differential  calculus  to  them.  Consider  ordinal  variables  x,  y,  and  z,  and  suppose  we  find  that  z 
varies  positively  monotonically  with  x  when  y  is  oonstant,  and  z  varies  negatively  monotonically 
with  y  when  x  is  constant.  Then  if  z  */tr,  y)  is  everywhere  differentiable  we  have  dzfbx  >  0  and 
dz/dy  <  0,  where  dz/3x  and  dz/dy  are  the  partial  derivatives  of  z  with  respect  to  x  and  y, 
respectively. 

In  general,  it  would  be  convenient  if  we  could  always  simplify  a  function  z~J{xj)  by  taking 
monotonic  transformations  of  our  variables  in  such  a  way  as  to  transform  the  function  to 
z*  *  x*xy*  or  z* »  x*  +  /,  where  z*.  x*,  and  /  are  the  transformed  variables.  Is  this  possible  in 
general?  It  is  not.  The  proof  is  not  difficult,  but  we  will  omit  it  here. 

But  the  unavailability  of  such  a  transformation  is  not  important,  since  if  we  restrict  ourselves 
to  qualitative  reasoning  and  ordinal  properties,  we  need  not  specify  any  but  the  ordinal  properties 
of  the  functions  we  deal  with;  the  exact  form  of  the  function  does  not  matter.  If,  in  the  last 
example,  x  undergoes  a  positive  monotonic  transformation  to  w,  so  that  now  z  *  g(w,y).  then  if 

4  >  0.  Sw  >  0. 


3.1 .  An  Example 

Consider  two  equations  connecting  a  pair  of  variables,  x  and  y,  with  an  exogenous 

parameter,  T,  in  the  first  equation; 
y  =A*T) 
fg(y) 

Suppose  the  system  described  by  these  equations  is  in  equilibrium  for  a  particular  value  of  T,  and 
that  T  is  now  displaced  by  an  amount,  iT.  We  now  solve  for  the  displacements  of  x  and  y  that 
restore  the  system  to  equilibrium,  dx  and  <ty  respectively; 

dy  *  fxdx  +  f-jdr 
dx*gydy 

where  the  subscripts  denote  partial  derivatives.  Then  we  find  that; 
frdT 


dy> 


U-4  gy) 
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dx-hhfL 


(1) 


Without  loss  of  generality,  we  can  take  fT  >  0,  for  we  can  always  replace  T  by  -T. 

IA:  Now  if  /,  and  gy  are  both  positive  and  fjgy  <  l.  then  dy/dT  >  0  and  dx/dT  >  0. 

IB:  lf/x  and  gy  are  both  positive,  but  fxgy  >  1,  then  dy/dT  <  0  and  dx/dT  <  0. 

HA:  Now  if/,  and  gy  are  both  negative  and  fg  <  l,  then  dy/dT  >  0  and  dx/dT  <  0. 

IIB:  If/,  and  gy  are  both  negative,  but  f£y  >  l,  then  dy/dT  <  0  and  dx/dT  >  0. 

Ill:  If/,  >  0  and  gy  <  0,  then  dy/dT  >  0  and  dx/dT  <  0. 

IV:  If/,  <  0  and  gy  >  0,  then  dy/dT  >  0  and  dx/dT  >  0. 

Table  3*1  summarizes  these  results.  Since  the  result  in  cases  IA  and  IB  depends  on  the 


cardinal  magnitude  of  fgy,  we  might  suppose  that  it  is  not  invariant  under  monotonic 
transformations  of  the  variables.  But  this  is  not  the  case.  Suppose  we  transform  x  monotonically 


to  w  *  p(x),  and  y  to  z  *  q(y).  Then 


and  = 

*  p  1  q 

y*  *7 


so  that 


_  Qyf iPx&y  _ 


Px<iy 


~fz8y 


The  product  of  the  partial  derivatives  of  the  two  dependent  variables  is  invariant  under  monotonic 


transformations  of  those  variables,  hence  the  sign  of  (1  -/,  gj  is  an  ordinal  property  of  the  system. 

Notice  that  the  conclusions  we  reach  by  solving  the  simultaneous  equations  are  not  always 
the  same  as  those  we  would  reach  simply  by  propagating  dT  through  the  equations,  step  by  step. 


For  consider  the  cases  examined  above: 

•  IA:  The  increase  in  T  causes  an  increase  in  y,  which  in  turn  causes  an  increase  in  x. 
If  we  trace  the  repercussions  further,  we  see  that  the  increase  in  x  will  cause  a 
further  increase  in  y,  so  that  there  will  surely  be  a  net  increase  in  x  and  y,  in 
agreement  with  the  algebraic  analysis  in  Case  IA.  Our  curiosity  may  be  aroused  as 
to  whether  the  process  will  converge  to  an  equilibrium  or  diverge,  a  question  that  did 
not  arise  explicitly  in  our  algebraic  treatment.  We  will  return  to  this  question  later.  IB 
is  similar. 

•  IIA:  The  increase  in  T  causes  an  increase  in  y,  which  causes  a  decrease  in  x.  The 
decrease  in  x  will  cause  a  further  increase  in  y,  and  so  on.  There  will  surely  be  a  net 
increase  in  y  and  a  net  decrease  in  x,  as  in  algebraic  Case  IIA.  But  again,  we  do  not 
know  whether  the  process  converges.  IIB  is  similar. 

•  III.  The  increase  in  r  causes  an  increase  in  y,  which  causes  a  decrease  in  x.  This 
time  a  decrease  in  x  will  cause  a  decrease  in  y,  and  this  will  cause  an  increase  in  x. 
In  this  case,  propagating  the  change  around  the  feedback  loop  makes  the  net  effect 
indeterminate,  and  again  we  do  not  know  if  the  process  converges. 

•  IV.  The  increase  in  T  causes  an  increase  in  y,  which  causes  an  increase  in  x.  Now 
the  increase  in  x  causes  a  decrease  in  y,  and  a  consequent  decrease  in  x.  As  in 
Case  III,  the  final  results  are  indeterminate. 
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Gin  we  aay  anytting  more  about  convergence?  Whan  we  propagate  tha  initial  impulse 
mound  tha  loop,  wa  have  -  ij&Ji  and  Ay*,  -  4 a x,,  ao  that  Ay*,  -  //, Ay.,  or 
Ay*^  «  (tyy*Ay^  Hanca,  tha  «um  £y,  wiM  converge  if  I//,!  <  l,  and  will  diverge  otherwise. 
By  introducing  an  explicit  dynamic  process  of  equiltoration  into  our  aaaumptiona,  wa  hava  arrived 
at  conditions  for  convarganca  of  tha  procaaa,  or  what  is  equivalent,  conditions  for  stability  of 
squiUbrium. 

Examining  Tabia  1,  wa  saa  that  in  tha  eases  of  convarganca  (tha  right-hand  half  of  tha 
table),  tha  net  affects  of  the  change  in  T  are  tha  same  in  sign  as  tha  gross  affects  of  tha  initial 
propagation  to  y,  and  from  y  to  x  Tha  net  change  in  y  will  always  hava  tha  same  sign  as  tha 
change  in  T,  while  tha  change  in  x  will  hava  tha  same  or  opposite  sign  as  gy  is  positive  or 
negative,  respectively.  Tha  sign  of  tha  feedback  link,  fx  is  irrelevant  -  at  most  it  magnifies  or 
damps  the  effect  of  the  initial  impulse  as  fxgy  is  positive  or  negative,  respectively. 


Net  Effects  of  Change  In  dT 

B 

A 

»*0y 

>1 

fiOy 

<1 

f* 

Oy 

dx 

dy 

dx 

dy 

l 

+ 

+ 

• 

- 

+ 

+ 

It 

- 

• 

+ 

- 

- 

♦ 

III 

♦ 

- 

• 

♦ 

- 

♦ 

IV 

- 

♦ 

+ 

♦ 

♦ 

♦ 

Table  3*1:  Propagation  Effects 


3.2.  Dynamics  and  Stability 

From  tha  results  of  tha  last  section,  wa  laam  that  even  in  vary  simple  systems  with 
fsedback  loops,  fits  conclusions  wa  can  draw  about  shifts  in  equilibrium  may  depend  critically 
upon  assumptions  of  system  stability.  But  the  conditions  for  stability  are  not  uniqueiy  determined 
by  tha  equilferium  aquations  alone.  They  depend  also  on  tha  process  of  adjustment  whan 
equilibrium  is  disturbed. 

Consider  again  tha  system: 

yJW) 

*•&) 

‘rh*  informal  argument  made  earlier,  that  fygt  <  l  was  a  necessary  and  sufficient  condition  for 
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stability,  was  derived  implicitly  by  embedding  our  equilibrium  equations  in  the  following  dynamic 
structure: 

iJ,) 

for  which  the  previous  pair  of  equations  is  the  (stable  or  unstable)  equilibrium.  Since  we  hold  T 
constant  after  the  initial  impulse,  dT,  we  may  simply  disregard  this  variable  and  consider  the 
behavior  of  the  system  moving  from  a  non-equilibrium  initial  position,  and  y0,  with  T  fixed.  We 
obtain  the  second-order  difference  equation: 

y,*2  -/W) 

For  small  disturbances  from  equilibrium,  we  can  replace  this  by  the  linear  system: 

y,+imAy, 

where  A~fx  gy.  Let 
so  that 

(AJty*2MA(Ajty, 

A*2*  1 

o- 

This  expression  for  yt  implies  that  y  will  converge  iff  A  <  1.  which  was  also  our  conclusion  in 
Section  3.1 .  In  general,  assuming  different  processes  of  adjustment  may  lead  to  different  stability 
conditions,  and  as  a  consequence,  to  different  net  effects  of  a  shift  in  parameters. 

3.3.  Propagation  and  Stability 

The  description  of  value  propagation  in  Section  3.1  is  very  similar  to  the  procedure  for 

solving  confluence  equations  to  derive  the  dynamic  behavior  of  a  system  employed  by  de  Kleer 

and  Brown  in  ENVISION  [1].  Consider  again  the  system  in  the  example. 
yfcT) 
x-«(y) 

If  we  assume  that  these  equations  represent  components  or  oonduits  that  transmit  signals 
instantaneously  (in  other  words,  that  the  equilibrium  relations  represented  by  these  equations  are 
restored  instantly  after  a  disturbance,)  we  can  derive  a  description  of  the  dynamic  behavior  of  the 
system  by  differentiating  the  equations  as  follows; 

£-/,§*/r  f  <*> 


« 
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<n 


As  in  the  example  in  Section  3.1 ,  we  assume  ^  >  0  and/r  >  0. 

Case  I:  Assume  /_.  $  >  0.  Equations  (2)  and  (3)  become  the  following  confluence 


equations; 


dy  *  dx  +dT 

(5) 

dx*dy, 

where  d»  denotes  the  sign  of  Given  the  input  disturbance, 

dT*+,  (6) 

the  values  of  dy  and  dx  can  be  determined  by  propagating  the  disturbance  through  the  confluence 
equations.  In  this  case,  we  have  a  choice  of  propagating  it  through  confluence  (4)  to  dx  or  to  dy. 
Case  IA:  Assuming  that  dx  is  initially  negligible,  we  propagate  the  input  to  dy  through  (4)  to 


obtain 


dy«+, 

Propagating  this  to  dx  through  (5),  we  have 
dx*+. 


Propagating  (8)  back  to  (7)  does  not  produce  a  contradiction.  Therefore,  (7)  and  (8)  are  a 
solution.  Furthermore,  if  we  now  assume  that  dT  is  negligible  and  propagate  (8)  to  dy  through  (7), 
the  value  obtained  for  dy  would  be  the  same,  +.  Therefore,  the  feedback  must  be  positive. 

Case  IB;  Likewise,  if  we  assume  that  dy  is  initially  negligible,  we  obtain 
dx  *  -  and  dy  *  -. 

Other  cases  lead  to  the  same  solutions,  as  shown  in  Table  3-1 .  Notice  that  in  solving  the 
problem  by  sign  propagation,  as  is  done  in  ENVISION,  it  is  the  assumptions  about  the  orders  of 
propagation  that  lead  to  different  solutions  in  Case  IA  and  IB  above.  On  the  other  hand,  in 
Section  3.1,  assumptions  about  the  value  of/xgy  lead  to  the  two  different  solutions.  Let’s  explore 
these  assumptions  in  terms  of  stability.  Equations  (2)  and  (3)  can  be  rewritten  as  follows: 

w 

Vrf  CO) 


f*8y<  1 


dy/dt,  dxjdt  >  0,  since  (l-/^p  >  0 ,fT  >  0,  and  dTIdt  >  0.  In  terms  of 
propagation  this  means  that  \  dy/dt  \  >  \fx(dx/dt)\,  hence  the  effect  of  dT/dt  on 
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dy/dt  dominates  the  effect  on  dx/dt.  This  is  equivalent  to  propagating  87  to  8y 
first,  i.e.  Case  IA  above. 

fay  »  1  dy/dt,  dx/dt  <  0.  In  terms  of  propagation  this  means  that  \dy/dt\  <  \fx(dx/dt)\, 

hence  the  effect  of  dT/dt  on  dx/dt  dominates  the  effect  on  dy/dt.  This  is 
equivalent  to  propagating  87  to  8z  first,  i.e.  Case  IB  above. 

As  shown  in  the  previous  section,  A  is  the  stable  system,  and  B  is  the  unstable  system, 
though  the  question  of  stability  cannot  be  answered  by  solving  confluence  equations.  Notice  that 
although  the  stability  assumption  is  weaker  than  the  assertion  that  8*  or  3y  is  negligible, 
nevertheless  it  is  strong  enough  to  determine  which  term  in  Equation  (2)  dominates. 

4.  Non-linear  Dynamic  Systems  with  Ordinal  Variables 

Thus  far  we  have  been  concerned  with  static  functional  relations  among  variables.  Now  we 
turn  to  time-dependent  relations  expressible  as  differential  equations.  The  discussion  will  be 
limited  to  the  simplest  case  of  first-order  non-linear  differential  equations  in  three  variables, 
including  time.  Any  such  system  can  be  written  as: 


(11) 

(12) 

Taking  the  ratio  of  Equations  (12)  to  (1 1),  we  also  have: 

4r  _«(*.>) 
d*  A*o>) 

(13) 

Here  we  are  on  more  familiar  ground,  for  textbooks  on  differential  equations  do  often 

present  the  kind  of  qualitative  treatment  of  non-linear  systems  that  we  are  about  to  describe  here. 

Hence,  we  can  be  brief.  In  the  x-y  plane,  we  can  draw  the  direction  field  that  describes  the  paths 

of  the  system  from  any  initial  conditions.  At  any  point  in  such  a  path,  the  slope  of  the  path  will  be 

given  by  equation  (13).  Figure  4-1  depicts  the  direction  field  for  (13)  in  the  case  where 
A *0»)  •-x  +  b  log(y)  and 

+  alog(x),  (14) 

If  we  set  j»0,  then  we  obtain  the  curve,  My)*0.  along  which  the  intersecting  paths  of  the 

direction  field  are  vertical.  Similarly,  if  we  set  %  *  0,  we  obtain  the  curve,  g(xjy 0,  along  which  the 

dt 

intersecting  paths  of  the  direction  field  are  horizontal.  The  points  where  —  *  £  =  o  are  the  (stable 

dt  dt 

and  unstable)  equilibrium  points  of  the  system.  In  Figure  4-1 ,  at  A,  where  the  paths  for  increasing 
t  converge  on  an  equilibrium  point,  the  equilibrium  is  stable;  at  B,  where  they  diverge  from  one,  it 
is  unstable.  At  saddle  points,  some  paths  may  converge  while  others  diverge,  making  the 
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Fig urt  4*1 :  Direction  Field  for  Equations  ( 1 4) 
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equilibrium  unstable. 

Sacks  developed  a  program  called  PLR  (Piecewise  Linear  Reasoner)  which  examines  the 
properties  of  the  direction  field  of  a  system  of  non-linear  differential  equations  to  determine  its 
qualitative  behavior  [7].  In  particular,  he  divides  the  space  into  regions  in  each  of  which  the 
system  can  be  approximated  by  a  linear  system.  In  each  such  region,  he  examines  the  direction 
field  to  determine  the  path  of  the  system  within  and  out  of  the  region.  The  global  behavior  of  the 
system  can  be  characterized  qualitatively  in  terms  of  the  system's  path  through  these  regions. 

Weld  has  proposed  a  technique,  which  he  calls  Comparative  Analysis,  and  which  can  draw 
some  conclusions  about  the  displacement  of  the  paths  and  of  equilibria  with  changes  in  system 
parameters  [9].  In  Section  6  we  will  show  how  these  sorts  of  conclusions  can  be  drawn  directly 
by  our  methods  of  ordinal  reasoning. 

In  non-linear  systems  of  differential  equations,  the  paths  of  the  direction  field  may  also  form 
closed  curves  in  the  x-y  plane,  so  that  the  system  describes  a  periodic  motion  around  the  curve. 
Such  paths  are  called  limit  cycles.  They  are  stable  if  neighboring  paths  all  converge  to  them 
(from  both  inside  the  closed  curve  and  outside),  unstable  otherwise. 

It  is  obvious  that  these  properties  -  the  distribution  of  equilibrium  points  and  limit  cycles,  of 
a  system  of  differential  equations,  and  the  stability  or  instability  of  these  equilibria  and  cycles  -- 
are  invariant  under  ordinal  transformations  of  x  and  y.  For  consider  such  a  transformation, 
x*  *  <Kx),  y*  =  v0>).  Since  the  orderings  of  x-coordinates  and  of  y-coordinates  are  not  altered,  all 
direction  paths  will  remain  intact.  If  they  converged  on  a  point  or  limit  cycle  in  the  x-y  plane,  they 
will  converge  in  the  x‘-y*  plane. 

Consider  the  case  where  the  system  defined  by  equations  (11)  and  (12)  has  an  equilibrium 


point  which,  by  transformations  of  the  coordinates  and  without  loss  of  generality,  we  can  place  at 
the  origin.  If  we  expand  /(•)  and  $(•)  in  Taylor's  series  and  neglect  the  higher-order  terms,  the 
system  is  approximated  by 

(15) 

dx  , 

—  •ax*  by 
dt 

(16) 

^-•cx*dy 

dt 

Integrating  these  equations,  we  obtain: 
x  *  A  e^M+ZJ 

(17) 
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ymCeW+Defa) 


where  X»  (o+d)±[(a+dj1-4(cd-bc)]lf2.  The  equilibrium  will  be  stable  iff  both  solutions  for  X  have 
negative  real  parts.  But  this  condition  will  hold,  in  turn,  iff  (a+d)  <  0  and  ( bc-od)  <  0.  >From 
Equation  (15)  (dy/dx)^^  ■  -ajb,  while  from  (16)  (dy/dx)ttyld^0*-c/d.  So  stability  depends  on 
which  of  these  two  slopes  is  the  greater.  If,  for  example,  a  and  d  are  negative,  while  b  and  c  are 
positive,  stability  requires  that  be  <  ad,  whence  -c/d  <  -a/b.  Then  the  slope  of  the  curve  for  — 

di 

must  be  steeper  than  the  slope  of  the  curve  for 

dt 


For,  by  the  chain  rule  for  differentiation, 

Subtracting  these  two  quantities,  we  have: 

<£>£.o  -  <£». 

dx  dt  dx  dt  rty  ^  dx  &  dx  & 

Since  the  first  two  factors  on  the  right  side  are  always  positive,  the  sign  of  the  difference  of 
the  slopes  of  the  transformed  variables  is  the  same  as  the  sign  of  the  difference  of  the  slopes  of 
the  original  variables.  Hence,  stability  is  preserved  under  positive  monotonic  transformations  of 

the  variables.  This  can  be  generalized  as  follows: 

Theorem  3:  For  a  system  of  differential  equations: 
dxi 

*/■  ■ '  ■  i^j)i  i*  li  ■  ■  i  a, 

the  stability  of  the  equilibrium  points  is  invariant  over  positive  monotonic 
transformations  such  as  wi  •  gt  (x,). 

Proof: 

The  eigenvalues  of  the  transformed  equations  are  the  same  as  those  of  the  original 
differential  equations. 

dwt  dwi  dxt 

~dr=dTi~diXg‘^ 

Linearizing,/;,  g?  about  the  equilibrium  point,  using  a  Taylor  series  expansion, 

/**  j 


“*«>)  +  •  • 

where  Iq*  (x10.  Xj0,  •  •  • ,  x^  at  equilibrium.  At  equilibrium/;  *  0,  and  with  no  loss  of  generality  let 
us  assume  x ^  0.  Therefore,  neglecting  higher  order  terms, 
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dw 

~d 

where 


dw-  ,  *  df- 

r«.  xy  -f/1  (wp  and  «,  -  (—\e. 


Substitute  for  xy  in  terms  of  wy  by  using  a  linearized  relation 


about  the  equilibrium  point. 

+  rVw'«‘H'»)  +  ••• 

Assume  wA  ■  0,  with  no  loss  of  generality.  Therefore, 

^  A 
-2TMZavwj 

J^dw  dx. 

where,  a,.  *  (— I-ri).  Notice  a,.  *  1. 

v  axt  aWj 

In  order  to  determine  the  eigenvalues  of  this  equation  consider  the  following  determinant. 


say  A. 

/i,-X,  a,2/12.  •  •  ■  0\j*-\f\A 
“21/21*  ‘  ■ 

=  0 

an\f»V  a«2/(i2'  '  *  •  In*  “*■ 

The  expansion  of  this  determinant  is  given  by 
1^1  -  X  a3T  '  -  a„v)IM 

a 

where  a  is  all  the  permutations  of  (1,  2,  ....  n),  and  |P0|»  ±  1  is  the  permutation  matrix  of  A. 
Therefore. 

'  ■  '  aJ  "Aafty&y  ' '  -4v 

where  */u  -  X  for  i  ■  8.  Notice  that  all  the  terms  with  non-repeating  indices  always  form  a 
cycle  with  each  index  in  both  positions,  i.e.  a,  _  and  an ..  Therefore  —  is  always  canceled  by 

This  implies  that  the  characteristic  polynomial  of  the  transformed  equation  is  the  same  as 
that  of  the  original  equation,  therefore  that  the  eigenvalues  are  exactly  the  same.  •  •  •  • 


5.  Extreme  Values 

Let  y«Ax)  be  a  single-valued  function  of  x,  and  x*  =  «x)  and  y'  =  v(>)  be  positive 
monotonic  transformations  of  x  and  y,  respectively.  Let  f[x)  have  a  local  maximum  at  xg,  so  that 
for  x  in  some  interval  about  x0,y(xrf*y0>y(x).  Now  y*>  y* iff  yi  >  >(x).  Hence,  /(xg)  >  y*(x*), 
where  xj  ■  9(xq)  and  x* »  $(x)  for  any  other  x  in  the  interval.  It  follows  that  maxima  of  the  original 
function  correspond  to  maxima  of  the  transformed  function. 

This  result  may  seem  counterintuitive,  since  whether  a  stationary  value  of  a  differentiable 
function  is  a  maximum  or  a  minimum  depends  on  the  sign  of  the  second  derivative,  and  this  sign 
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is  not.  in  general,  invariant  under  positive  monotonic  transformations  of  the  variables  (concavities 
can  change  to  convexities,  and  vice  versa).  However,  it  is  easy  to  show  that  the  sign  of  the 
second  derivative  is  invariant  in  the  neighborhood  of  a  stationary  point.  The  proof,  which  again 
makes  use  of  the  chain  rule  for  differentiation,  is  straightforward,  and  will  be  omitted  here. 

6.  Comparative  Statics 

Much  qualitative  analysis  takes  the  form  of  comparative  statics.  The  (stable)  equilibrium  of 
a  system  is  displaced  by  a  change  in  one  of  its  parameters,  and  we  wish  to  know  how  the  values 
of  the  system  variables  are  changed  when  it  settles  to  its  new  equilibrium.  In  many  cases  of 
interest,  the  shapes  of  the  system  functions  are  not  known,  but  only  the  signs  of  their  derivatives; 
and  only  the  sign  of  the  disturbance,  and  not  its  magnitude,  is  given. 

To  see  what  kinds  of  conclusions  can  be  reached  under  such  circumstances,  let  us 
consider  a  simple  example  from  economics.  (See  Figure  6-1.)  For  a  certain  commodity,  the 
quantity,  q,  will  be  supplied  by  producers  if  the  price  is  ps  -  p^q),  and  the  same  quantity  will  be 
purchased  by  consumers  if  the  price  is  pD  *  pD(q).  More  will  be  supplied  if  the  price  is  higher 
(~  >  0),  and  less  will  be  purchased  if  the  price  is  higher  >  0).  The  market  will  be  in 

4q  dq 

equilibrium  when  p^qrf  ■  pD{q$  -  Po- 

Now  suppose  that  a  sales  tax  is  imposed  on  the  commodity,  so  that  the  r.et  supply  price  for 
any  quantity  is  increased  by  the  amount  of  the  tax.  Call  the  new  supply  price  ps  *  p^q)  >  p^q), 
and  p^o)  *  PD(?o)  *  Pa- 

Theorem  4:  The  equilibrium  price  will  be  increased,  p0  >  p0,  and  the  equilibrium 
quantity  will  be  decreased,  ?0  <  q#  Moreover,  since  the  tax  is  such  that  p^q)  increases 
with  an  increase  in  q,  the  equilibrium  increase  in  price  will  be  smaller  than  it  would  have 
been  if  the  quantity  exchanged  had  remained  constant: 

P^A o)  mPo 

Proof:  These  results  do  not  depend  on  continuity  or  differentiability  of  the  functions.  For 
all  q  >  qQ,pD(q )  <  Pq,  since  the  demand  price  decreases  when  the  quantity  increases.  But 
p£q)  <  pQ,  since  the  supply  price  increases  under  the  same  conditions. 

For  all  q,  p£q)  >  p^q),  the  difference  being  the  amount  of  the  tax.  Therefore, 
PD(q)  <  Po  <  P^4)  <  Pfa)  fof  all  q  >  <7o>  and  pD(q)  <  p^q)  in  this  range  of  q.  Hence,  if  a  new 
equilibrium  exists,  we  must  have  p0  <  qQ.  Now  for  all  q  <  q0  (including  5q),  p^q)  <  p^q^},  hence 
p&o)  >  Ay(?o)  mPo-  Similarly,  pD(q)  >  pD(qJ,  so  that  pD(fy  *  P0  >  Po  ■ 
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7.  Conclusion 

In  these  pages  we  have  sketched  out  some  of  the  properties  of  functions,  and  especially 
continuous  differentiable  functions,  that  are  invariant  under  monotonic  transformations  of  the 
variables,  and  have  shown  how  these  properties  can  be  used  to  analyze  phenomena  where  the 
variables  are  only  defined  ordinally. 

The  main  idea  is  that,  in  the  case  where  two  variables  are  positively  (negatively) 
•  monotonically  related,  their  positive  monotonic  transforms  are  positively  (negatively) 

monotonically  related.  If  x  is  a  monotonic  function  of  y  and  y  of  z,  then  x  is  a  monotonic  function  of 
z.  if  we  assign  the  number  zero  to  a  monotonic  relation  if  it  is  positive,  and  1  if  it  is  negative,  and 
if  xpZj.x^xj,  ...  x.  are  a  sequence  of  variables,  each  of  which  is  a  monotonic  function  of  its 
successor,  than  the  sign  of  the  relation  between  x,  and  x„  will  be  positive  or  negative  as  the  sum, 
modulo  2,  of  the  numbers  assigned  to  the  intervening  relations  is  0  or  1 . 

Much  reasoning  about  magnitudes,  even  in  the  sciences,  takes  place  without  any  help  from 
mathematical  formalisms,  but  solely  in  terms  of  ordinary  language.  As  the  illustrations  in  this 
paper  suggest,  a  good  deal  of  this  reasoning  makes  implicit  use  of  the  properties  of  ordinal 
variables  and  monotonic  transformations.  In  both  the  social  science  and  natural  science 
literatures,  much  reasoning  is  also  done  from  diagrams,  without  concern  for  the  exact  forms  of 
the  functions  depicted  or  the  cardinal  values  of  variables.  This  reasoning  also  makes  use  only  of 
the  ordinal  properties  of  the  magnitudes  depicted  in  the  diagrams.  We  would  conjecture  that 
much  if  not  most  everyday,  'commonsense*  reasoning,  when  it  is  concerned  with  magnitudes,  is 
-  actually  implicit  reasoning  about  ordinal  relations. 
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